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Notice 


This report was prepared for the Building and Fire Research Laboratory 
of the National Institute of Standards and Technology under grant number 
60NANB5D0113. The statement and conclusions contained in this report 
are those of the authors and do not necessarily reflect the views of the 
National Institute of Standards and Technology or the Building and Fire 
Research Laboratory. 


Abstract 


A new type of near infrared (NIR) fire detector which utilizes a statistical analysis of apparent 
source temperatures of fires was developed. The NIR fire detector measures the spectral radiation 
intensities emanating from fires at 900 and 1000 nm. These measurements are used to obtain and a 
time series of apparent source temperatures. The near infrared radiation characteristics of five standard 
test fires specified m the guidelmes of the European Committee for Standardization were measured and 
utilized to develop a fire detection algorithm. The five test fires specified m the guidelmes imvolved a 
heptane pool, a polyurethane foam, a wooden crib, a smoldermg cotton and a smoldermg wood. The 
normalized power spectral density of the spectral radiation mtensity at 900 or 1000 nm combmed with 
the probability density function of the apparent source temperatures are sufficient to determme the 
presence of the test fires m the vicmity of the detector. The NIR fire detector was evaluated for 
sensitivity to both direct and reflected radiation from the five test fires. When direct radiation is 
incident on the NIR fire detector, four of the five test fires were always detected. The smoldermg 
wood fire was detected only if the NIR fire detector was very close. The NIR fire detector can also 
detect the heptane pool, the polyurethane foam and the wooden crib fires from reflected radiation. The 
response time of the NIR fire detector ranged from approximately 1 to 4 mmutes for the open fires. 
The NIR fire detector had no instances of false alarms during the test period. The response time can be 
shortened considerably (to a few seconds) by usmg a programmable DSP chip rather than a personal 
computer. We continue to seek SBIR fimding for developimg this faster version. 


A numerical technique based on the Discrete Probability Function (DPF) method m 
conjunction with a ray tracing algorithm was developed to handle radiative heat transfer problems 
im rectangular and cylmdrical enclosures without participatmg media. The numerical technique 
has potential to assist buildig designers to find the optimal detector location for multiple room 
monitormg. The DPF method is much more accurate than conventional Monte Carlo methods for 
solving the surface radiative heat transfer problems studied im this project. The numerical 
evaluation confirmed that the detector can successfully detect fires from reflected radiation if its 
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Description 

Velocity of light, m/s 

Spectrum function, s 

Normalized spectrum function 
frequency 

Probability density function of variable o 
Planck’s constant 

Radiation intensity 

Boltzmann constant 

Length of the tube 

Samplmg frequency 

Number of bins m the DPF of position 
Discrete probability function of variable o 
Radius of the tube 

Auto-correlation coefficient 

Direction cosines 

Radial coordinate 

time 

Temperature, K 

Axial coordinate 


Azimuthal angle 

Cosme of the polar angle 

Polar angle 

Hemispherical reflectivity 
Hemispherical specular reflectivity 
Hemispherical diffuse reflectivity 
Wavelength 


blackbody 

subscripts for bin values 

Wavelength 

First measurement wavelength (1000 nm) 
Second measurement wavelength (900 nm) 
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1. Introduction 


The ability to detect fires m residential and industrial venues has been and will continue to 
be of the utmost importance. In 1995 alone, public fire departments responded to 1.965 million 
fires (Krater, 1996). 573,500 of these fires were structure fires resulting in $7.62 billion dollars in 
property damage. Seventy four percent, or 425,500, of these fires occurred in residential 
properties causmg $4.363 billion m damage. More important than monetary concems, fires also 
claimed 4,585 civilian lives and injured an additional 25,775. The number of fire deaths in the 
home rose 6.3% im 1995 to 3,640, accounting for 79% of all fire deaths. An additional concer is 
false alarms. It is estimated that 95% of all alarms from smoke detectors in the U. S. during the 
1980’s were unnecessary (Hall, 1989). 


New fire detection concepts and algorithms are justified only if they improve upon existing 
ones with lower false alarm rates and greater sensitivity to starting fires. In addition, the detectors 
and signal processmg mstruments should be easy to operate and maintain, have high flexibility and 
be relatively mexpensive (Luck, 1992). Currently residential fire detectors include optical smoke 
sensors, ionization smoke sensors and temperature sensors (Grosshandler, 1992). 


Conventional smoke sensors utilize light scattermg or smoke ionization measurements to 
detect a fire, while temperature sensors utilize thermocouple measurements. The are three 
disadvantages with conventional smgle sensor detectors: (1) there is a significant time delay 
between the start of the fire and the transport of the combustion products to the location where 
the detector is mounted; (2) m instances when there are impermeable barriers (such as smoldermg 
mside walls), the fire is not easily detected even m advanced stages, and; (3) smgle sensor 
detectors mvolve a high rate of false alarms due to changes m the operatmg environment. 
Combinations of smoke sensors and odor sensors which mvolve multiple fire signatures are less 
prone to false alarms (Okayama et al, 1994). However, multiple sensors mvolve greater 
construction cost and increased complexity of signal processing hardware and software. 


More recently, there has been increased interest im the use of radiation emission sensors 
(flame detectors) as an alternative to smoke and heat sensors (Middleton, 1989). The three major 
advantages of emission sensors are: (1) their ability to survey the entire room for fire mitiation, (2) 
their fast response time, and (3) false signals can be readily distinguished smce most fires are 
unsteady with unique frequency content, leading to unambiguous discrimmation based on the 
power spectral density of the measured mtensities (Grosshandler, 1992). 


Single channel flame detectors operate either im the ultraviolet (where solar radiation is 
totally absorbed by the earth’s atmosphere) or in the infrared (where flame emission is primarily 
from hot CO) parts of the spectrum. Ultraviolet signals from flames are normally very low 
leading to false alarms from indoor radiation sources such as incandescent lights, arc welding 
processes, etc. Therefore, ultraviolet sensors are limited to outdoor usage where interfermg solar 
radiation is absorbed by the earth’s atmosphere. Another disadvantage of ultraviolet flame 
detectors is that any contamination of the optical windows causes a significant loss of sensitivity. 


Infrared flame detectors are used for large indoor areas such as aircraft hangars and 
warehouses where direct solar radiation is mmimal. Single channel mfrared detectors look for 
radiation emitted from hot CO, gases present in most flames at wavelengths around 2.7 um or 4.4 
um. These single channel detectors have precision band-pass optical filters m front of them to 
detect fires while successfully rejecting solar radiation. The major problem with single channel 
detection is that since only one channel of information is present, the chances of false alarms are 
relatively high (Okayama et al., 1994; Middleton, 1989). 


The false alarm problems present with single channel detection can be partially alleviated 
by using two channels of information. Typically two-channel flame detectors use one channel m 
the infra-red (typically at 4.4 um) to detect hot combustion products. The second channel is 
chosen above or below the 4.4 um band where there is a high level of solar radiation coupled with 
low levels of flame radiation. The addition of the second channel is purely for the prevention of 
false alarms by rejecting interference (such as direct solar radiation) from a continuum source that 
does not have the ubiquitous 4.4 um CO band. Fire is still detected usmg the 4.4 um mfrared 
channel, and in cases where fire is present along with the mterfermg source, it might be difficult to 
resolve the signal unambiguously (Middleton, 1989). Further, highly lummous fires may go 
undetected if the detector is tuned to non-lummous fires. Commercial production of smgle 
channel infrared flame sensors which are msensitive to solar radiation, or a combination ultra- 
violet/mfrared and even two-channel mfrared flame sensors has been initiated for use m imdustrial 
applications (Middleton, 1989). 


A fiber optic fire sensor that uses correlation between radiation at two wavelength spaced 
far apart m the visible to 2.0 micron wavelength band has been recently utilized to detect diesel 
fires (Wetzork et al, 1992). The fiber optic fire sensor utilizes the high degree of correlation 
between the intensities at the two wavelengths to provide immunity to false alarm sources. 
However details regardmg the hardware of the fiber optic fire sensor or the fire detection 
algorithm were not reported. 


The two most distinguishing features of a natural fire, particularly a luminous one, are its 
apparent source temperature and the power spectral density of the radiation intensities emitted 
from it. Two-wavelength pyrometric measurements conducted in luminous pool and jet fires 
(Sivathanu and Faeth, 1990; Sivathanu et al, 1991; Sivathanu and Gore, 1991) indicate that the 
peak temperatures within these fires are m the range 1400 + 300K. The power spectral density of 
natural fires show a wide range of frequency present in them. 


The additional advantage of utilizmg radiation for fire detection is that the sensor can be 
placed away from the direct view of the fire. This allows multiple room monitoring with a single 
fire detector. However in arbitrarily shaped residential and commercial buildings, it is difficult to 
determine the right location for the placement of the fire detector. Typically, the fire detector 
should be located in places where the intensity incident on the detector is at a maximum. Design 
procedures for optimal placement of fire detectors have not been developed so far. For 
evaluation and eventual utilization of infrared radiation based fire detectors in residential and 
commercial buildmgs, a design tool is needed. 


The calculation of radiative transfer in enclosures with and without participating media is 
of particular interest in developing such design tools for all radiation based fire detectors. 
Solutions to the equation of radiative transfer m enclosures have been obtained using a wide 
variety of methods. Pioneering computations (Hottel, 1954; Oppenheim, 1956; and Sparrow et 
al, 1961) mvolved solutions to the radiative transfer equations for diffuse surfaces. Eckert and 
Sparrow (1961) pioneered the use of the angle-factor method for analyzing radiative transfer in 
enclosures with specular surfaces. Subsequently different approximate methods based on 
extensions of the angle factor or the script-F method have been used to analyze radiative transfer 
in enclosures having specular and diffuse surfaces (Sparrow et al., 1962; Bobco, 1964; Bevans 
and Edwards, 1965; and Sarofim and Hottel, 1966). In addition, analytical methods for the 
calculation of radiative heat transfer for simple surfaces such as cylindrical and conical cavities 
(Lim and Sparrow, 1965) and axisymmetric passages (Rabl, 1977; and Mahan et al., 1979) have 
been developed. 


Statistical methods for the calculation of radiative heat transfer in enclosures using Monte 
Carlo simulations have been published a few decades ago (Howell and Perlmutter, 1964; and 
Perlmutter and Howell 1964). Initially, unbiased or Direct Monte Carlo (Corlett, 1966; and Toor 
and Viskanta, 1968) simulations were used to calculate the thermal radiation m enclosures with 
mixed specular/diffuse surfaces. A review of the earlier Monte Carlo simulations of radiative 
transfer problems is found m Howell (1968). 


Direct Monte Carlo methods are computer mtensive, and therefore fully vectorized or 
parallelized codes are very helpful in reducing the costs associated with long run times. Monte 
Carlo algorithms for an irregular enclosure with arbitrary values for surface properties 
(transmitivity, specular and diffuse reflectivity) have been implemented usmg fully vectorized 
codes (Bums and Pryor, 1989). The accuracy and convergence of direct Monte Carlo algorithms 
have been documented (Maltby and Burs, 1991). A grid shading algorithm (Burns and Pryor, 
1989) was shown to reduce computational time significantly. For complex 2-D and 3-D heat 
transfer problems, similar changes in the photon tracking kernel are necessary to obtaim reasonable 
computational efficiency, and the authors (Bums and Pryor, 1989) suggested a weighted moving 
average process to reduce the number of photon trajectories required m such calculations. 


There are several methods of reducing the sample size in Monte Carlo methods such as 
importance sampling and correlation and statistical estimation techniques (Kahn and Marshall, 
1953). For radiative transfer calculations, weighting and biasmg mdividual photons (Lanore, 
1971) and importance sampling in conjunction with angular discretization (Burgart and Stevens, 
1970) have been shown to reduce the computational time and imcrease accuracy for deep 
penetration problems. 


Other schemes for improving accuracy such as the energy partition method (Shamsundar 
et aL, 1972; and Modest, 1978) have been used to combine deterministic elements with a Monte 
Carlo simulation of radiative transfer calculations. Haji-Sheikh (1988) and Burns et al. (1990) 
provide a detailed literature review of the more recent Monte Carlo simulations of radiative 
transfer. Haji-Sheikh (1988) provides a very thorough exposition of the importance sampling 
method in radiative heat transfer problems for variance reduction. Monte Carlo methods are 


currently being used in a variety of problems demonstrating their applicability for irregular shapes 
(Parthasarathy et al, 1994), discrete arrays (Drost and Welty, 1992), diffuse and specular 
reflectance and transmittance in axisymmetric enclosures (Bums et al, 1992). Recent papers 
(Tong and Skocypec, 1992) have shown that application of similar Monte Carlo approaches to a 
bench-mark problem can result in differences m solutions by a factor of 2 to 3. 


Markov chain theory (Naraghi and Chung, 1984; and Billings et al, 1991) has been used 
to calculate the radiation interchange in enclosures with a finite number of isothermal diffuse and 
specular surfaces. Similar to the exodus modification of the Monte Carlo method (Emery and 
Carson, 1968), Markov chain theory does not need random number generators, elimmatmg one 
source of statistical fluctuations in the solutions. For smgle parameter state defimition of the 
Markov chains, the calculation of the transition probability matrix for multiply (more than twice) 
reflected photons or for arbitrarily varymg specular reflectivity proves to be challengmg. A two 
parameter state definition m conjunction with the Markov cham method (Billings et al, 1991) 
makes it easier to treat specular reflections m enclosures. Billmgs et al (1991) modeled the 
radiative heat transfer im a 2-dimensional square enclosure usmg this method. The biggest 
advantage of the Markov chain method is the ability to reuse the probability transition matrix. 


Recently, to address the disadvantages associated with the use of random number 
generators im the Monte Carlo methods, we developed a Discrete Probability Function (DPF) 
method to calculate the spectral radiation mtensity leaving a representative path m a flame 
(Sivathanu and Gore, 1993). The DPF method mvolves discretization of the probability density 
functions (PDFs) of the initial and boundary conditions of a variable mto a finite number of bins 
(histogram representation). The governmg equations are then applied to the individual bin values. 
Application of the governing equations to mdividual bi values is similar to their application to 
randomly selected values m the Monte Carlo method. In the DPF method, the values and 
probabilities associated with the solutions of the governmg equations obtained are rebimned to 
provide the PDF of the solution. In the Monte Carlo method, the random sampling algorithm 
selects more of the higher probability values and less of the lower probability values of the imitial 
and boundary conditions and source functions. The application of the governing equations to the 
resulting realizations leads to solutions that upon binning yield their PDF. 


The DPF method was subsequently applied to the calculation of the view factors in a 
cylindrical tube with diffusely reflectimg surfaces (Sivathanu and Gore, 1994). The computational 
cost of carrymg the mformation concerning the probabilities was found to be less than that of 
considermg a large number of random realizations required for Monte Carlo calculations in these 
two applications (Sivathanu and Gore, 1993, Sivathanu and Gore, 1994). 


Based on the background information provided above, the objectives of the present 
project were: (1) to characterize the near mfrared radiation intensities emanating from standard 
test fires, and (2) to develop numerical methods for the calculation of spectral radiation intensities 
incident on detectors from fires m arbitrarily shaped enclosures. The experiments and numerical 
simulations will be used to evaluate a near infrared fire detector, specifically suited for monitormg 
multiple room enclosures. 


2. Experimental Evaluation 


2 Principle of Operation 


The principle of operation of the near-infrared fire detector is similar to the two- 
wavelength optical pyrometer used for determmimg soot volume fractions and temperatures in 
laboratory scale fires (Sivathanu et al., 1991). The spectral radiation intensity emitted by a source 
at any wavelength can obtained from the equation of radiative transfer as: 


I, =e,hyp (2.1) 
where I, 1s the blackbody intensity at the unknown source temperature T and ¢, is its apparent 
spectral emissivity. 


For fire detection, the exact temperature and the emissivity of the source is not of great 
importance. Rather, the existence of high temperatures m the vicmity of the detector is sufficient 
to indicate the presence of a fire. Therefore, the spectral emissivity of any radiation source (either 
direct or reflected) can be assumed to vary mversely with wavelength irrespective of its chemical 
composition, or the spectral reflectivity of mtervenmg material Usmg this assumption, the 
apparent temperature of any source, determmed from the measured spectral radiation mtensities at 
two wavelengths can be defined as: 
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where h is the Planck’s constant, k is the Boltzmann constant, and c is the speed of light. The 
advantage of using two wavelengths close to one other is that the assumption of 1/A dependence 
for the emissivity of the source does not mtroduce a very large error on the apparent source 
temperatures obtained using Eq. (2.2). 


A schematic diagram of the near-infrared (NIR) fire detector is shown in Fig. 1. The fire 
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Figure 1. Schematic diagram of the NIR fire detector. 


detector consists of 90° view angle optics that collect and collimate the radiation intensity incident 
on it. The collimated radiation intensity is split into two parts which are mcident on two photo 
multiplier tubes (PMTs) that have narrow band pass filters (10 nm FWHM) centered at 900 and 
1000 nm im front of them. The output voltages of the two PMTs were monitored usmg an A/D 
board and a laboratory computer. These voltages were converted to spectral radiation imtensities 
using calibration constants obtained with a reference blackbody maimtamed at 1100 K. 


The peak spectral radiation mtensities in typical lummous fires such as acetylene and 
ethylene are at approximately 1500 nm (Gore, 1986). The mtensities at 1000 nm are 
approximately a factor of three to four lower than at 1500 nm However, there are three 
advantages to utilizmg the spectral radiation intensities at 900 and 1000 nm rather than m the 
infrared regions to obtain apparent source temperatures: (1) near infrared fire detectors are less 
expensive, more sensitive and less prone to degradation than infrared detectors, (2) the sensitivity 
of the two-wavelength pyrometric technique to temperature is more when the wavelengths are 
close together and well shifted from the peak of the Planck function, and (3) buildmg materials are 
much more reflective at 900 and 1000 nm than at 1500 nm, making detection of reflected 
radiation easier. The lower sensitivity of the mfrared detectors could be alleviated to some extent 
by utilizmg broader bandwidth optical filters in lieu of the 10 nm filter used here. However, the 
detector will have higher false alarm rates due to the lower rejection of environmental radiation. 
For the characterization studies reported m the followmg, the choice of 900 nm and 1000 nm 
wavelengths is close to the optimal, providing both a higher sensitivity to temperatures and an 
increased immunity to false alarms. 


The spectral radiation mtensities at 900 nm and 1000 nm incident on the two photo 
multiplier tubes were measured at 100 Hz. From the spectral radiation intensities, a time series of 
apparent source temperatures was obtained usmg Eq. (2). The time series of spectral radiation 
itensities and apparent source temperatures were analyzed to obtain their Probability Density 
Functions (PDF) and Power Spectral Densities (PSD). 


22 Description of Standard Test Fires 


Fire detectors commonly used m residential and commercial buildings are tested using six 
standard fires specified i the European Committee for Standardization (CEN, 1982) guidelines. 
Five of these standard fires are lummous. The sixth one is a non-luminous alcohol fire which 
cannot be detected by the present fire detector. The near-infrared radiation characteristics of the 
five lumimous fires were studied to develop an effective fire detection algorithm. 


A brief description of the five test fires is included in the following section. Some of the 
dimensions of the fires specified m the CEN (1982) guidelines were scaled down so as to be able 
to accommodate them in our laboratories. The open cellulosic fire (designated TF1 im the CEN 
(1982) guidelmes) consisted of seven stacks of beechwood sticks arranged to form a wooden crib. 
Each stack consisted of three beechwood sticks of dimension 1 cm x 2 cm x 25 cm The 
beechwood sticks were arranged in a cross layer pattern to form a square wooden crib 250 mm on 
each side. The wooden crib fire was ignited with 5 cm’ of methylated spirit placed at the center. 


The smoldering pyrolysis fire (designated TF2) consisted of 24 dried beechwood sticks 
arranged along the twelve radii of a grooved hot plate 220 mm im diameter. The beechwood 
sticks were 1 cm x 2 cm x 3.5 cm, and the hot plate was powered using a variable voltage 
controller so as to reach of temperature of 600 ° C in 11 minutes. The wood smolders at 
approximately 8 minutes. The radiation from the hot plate is very high and therefore a radiation 
shield was used to prevent the detector from directly viewmg the hot plate. During the 
characterization tests, which lasted approximately 3 mmute, a flame was not present. 


The glowimg smoldering fire (designated TF3) consisted of 90 cotton wicks, 80 cm long, 
fastened by a wire rmg 10 cm m diameter at one end. To maintain a cylindrical shape for the 
cotton wicks, they were wetted with water and allowed to dry m a current of air, while a small 
tension was applied at the other end. The cotton wicks were ignited using a propane torch. The 
flame were put out immediately after ignition, and the cotton wicks continued to glow for a long 
time. 


The open plastics fire (designated TF4) consisted of 3 mats of soft polyurethane foam 
stacked one on top of another. The mats were 15 cm im diameter and 25 mm m height. The 
polyurethane foam has not been treated with any flame retardmg additives and therefore is very 
flammable. The bottom mat was ignited at one comer and the flame spreads very rapidly, with 
the entire fuel bemg consumed m approximately 1 minute. 


The liquid fire (designated TF5) consisted of a n-heptane pool fire stabilized on a water- 
cooled 15 cm diameter stamless steel burner. The height of the stamless steel burner was 10 cm, 
and the flames were stabilized with a lip height of 10 mm. The heptane pool fire reached a steady 
state operating condition m approximately 3 to 4 mimutes. 


2.3. | Expermental Arrangement 


The near infrared fire detector was evaluated using both direct and reflected radiation. The 
location of the detector with respect to the fires for evaluation usmg direct and reflected radiation 
are shown in Figs. 2 and 3 respectively. For the experimental evaluation usig direct radiation, 
the standard test fires were placed at a distance of 2 meters from the NIR fire detector. 


Fire Fire detector 


Figure 3. Experimental arrangement used for evaluation using direct radiation. 
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Figure 3. Experimental arrangement used for evaluation using reflected radiation. 


For the experimental evaluation of the NIR fire detector utilizmg reflected radiation, the 
direct radiation was blocked using a shield. The NIR fire detector and the test fires formed the 
two vertices of an equilateral triangle as shown in Fig. 3. The radiation from the fire was reflected 
onto the NIR fire detector using a poster board. The effect of different materials was briefly 
investigated using a brick wall, and a window pane at the same location of the poster board. 
These results are also presented m the followimg. 


2.4. Results of the Experimental Evaluation 


The results of the experimental evaluation are provided separately for the three open fires 
(TF1, TF4 and TF5) and the two smoldering fires (TF2 and TF3). The last fire (TF6) is an non- 
lummous alcohol fire and cannot be detected by the present NIR fire detector. 


The probability density functions (PDFs) of spectral radiation mtensity at 900 nm 
emanating from the open fires as measured by the NIR mfrared detector are shown m Fig. 4. The 
PDFs of spectral radiation mtensity mcident on the fire detector when the fires were im its direct 
and indirect view are shown im the top and bottom panels of Fig. 4. When the open fires were m 
the direct view of the detector, the spectral radiation mtensities incident on the detector ranged 
from a mean value of 100 W/m’-pm-str. for the polyurethane foam to 300 W/m’-um-str. for the 
heptane pool When the fire detector is not m the direct view of the fires, the only radiation 
icident on it was that reflected by the poster board. The poster board has very low reflectivity. 
Therefore, the spectral radiation mtensities mcident on the fire detector were approximately 40 to 
50 times lower than those obtained from direct radiation. In addition, the open polyurethane foam 
and wooden crib fires were transient m nature. This causes the actual intensities measured by the 
detector to vary depending on the time lag before the measurement and the initiation of the fire. 
Therefore, the fire detection algorithm used im radiation detectors has to account for this widely 
varying magnitude in actual radiation intensities incident on it, without causing any false alarm. 
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Figure 4. PDFs of spectral radiation intensities from the open fires. 


The PDFs of spectral radiation imtensity at 900 nm emanating from the smoldermg fires 
are shown in Fig. 5. When the smoldering cotton fire is m the direct view of the detector, the 
spectral radiation imtensities incident on the detector ranged from a 0.25 to 0.5 W/m’-um-str. The 
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Figure 5. PDFs of spectral radiation mtensities from the smoldermg fires. 


spectral radiation intensities from the smoldermg wood fires were even lower, rangmg from 0.1 to 
0.18 W/m’-pm-str. Similar to the open fires, the intensities obtained from reflected radiation were 
much lower than those obtamed from direct radiation for the smoldermg fires. The spectral 
radiation intensities from reflected radiation for the smoldermg wood fire were below the 
detection limit of the NIR fire detector. 


The major observation to be noted from the PDFs of spectral radiation mtensities shown m 
Figs. 4 and 5 is that the magnitude of mtensities mcident on the detector can vary by 
approximately 5 orders of magnitude. Therefore, a fire detection algorithm based on the 
magnitude of mtensities will either miss some of the fires (if the sensitivity of the detector is low) 
or suffer from false alarms from environmental radiation (for high sensitivity detectors). 


The probability density functions (PDFs) of apparent source temperatures estimated from 
measurements of direct and reflected spectral radiation mtensities at 900 and 1000 nm are shown 
in Fig. 6. The PDFs of apparent source temperatures obtained with the three open fires (heptane, 
foam and wood) m the direct view of the detector are shown m the bottom panel of Fig. 6. The 
apparent source temperatures for the three open fires ranged from 1300 to 1800 K. 
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Figure 6. The PDFs of apparent source temperatures for the open fires. 


The PDFs of apparent source temperatures for the three open fires when only reflected 
radiation was incident on the detector are shown im the top panel of Fig. 6. For liquid pool fires, 
the highest temperatures are usually associated with very low soot volume fractions, resulting in 
comparatively lower spectral radiation mtensities. In addition, at the lowest temperatures, the 
spectral radiation mtensities are very low due to the exponential nature of the Planck Function. 
Therefore, after one reflection, only the middle range of temperatures (with some spectral biasing) 
are detected by the NIR fire detector. 
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For the wooden crib and polyurethane foam fires, mformation on the correlation between 
local intensities and temperatures is not available. Therefore, the changes m the shape of the 
apparent source temperature PDFs could be due to a combmation of factors, mcluding spectral 
biasing of the reflected intensities, shape of the flame, and correlation between local temperatures 
and emissivities. In addition, these fire were transient m nature, and the temperatures obtamed 
from direct and reflected radiation could be at different stages m the development of the fire. 
Despite the differences in the PDFs obtained from the direct and mdirect viewing of the fire, all of 
the temperatures fall withm 1000 to 2000 K. 


The PDFs of apparent source temperatures estimated by the NIR fire detector for the 
smoldering fires (cotton fiber and wooden pellets) are shown in Fig. 7. Smolderig fires have 
much lower temperatures than open fires. The PDFs of apparent source temperatures obtained 
from direct radiation for the smoldering cotton fire varied from 900 to 1100 K, and that for the 
smoldering wood fire varied from 700 to900 K, as shown m the bottom panel of Fig. 7. 


The PDF of apparent source temperatures estimated by the NIR fire detector usmg 
reflected radiation from the smoldermg cotton fire was approximately 200 K lower. In addition, 
the PDF was not continuous. This is because the reflected mtensities incident on the NIR fire 
detector from the smoldering cotton fire were very low. Therefore, the discretization error of the 
A/D converter became a significant factor m the apparent source temperatures estimated by the 
NIR fire detector. For the smoldermg wood fire, the reflected intensities were below the 
detection threshold of the NIR fire detector. 


Based on the PDFs of apparent source temperatures estimated by the NIR fire detector 
from the open and smoldermg fires, one criteria for the existence of the fire m the vicmity of the 
detector can be identified. This criteria is that at least 40% of the apparent source temperatures 
estimated by the NIR fire detector should be between 600 and 2500 K. This is admittedly a very 
broad range of apparent source temperatures. However, there are very few sources of 
environmental radiation that would have an apparent source temperature that is within this range. 
For mstance, solar radiation would yield an apparent source temperature of 6000 K. Two 
common source that could have apparent source temperatures between 600 and 2500 K are 
commercial hot plates and mcandescent bulbs. Therefore, usmg only the PDF of apparent source 
temperatures as a fire detection algorithm would result m false alarms m the presence of a hot 
plate or an mcandescent bulb. 


The second characteristics of natural fires that can be used to eliminate these sources of 
false alarms is the power spectral density of spectral radiation intensities. The power spectral 
density of spectral radiation mtensities emanating from natural fires have a very well defined low 
frequency component, and this can be exploited m fire detection algorithms. 


The normalized power spectral densities (PSDs) of spectral radiation mtensities at 1000 


nm for the three open fires are shown m Fig. 8. For a unsteady signal I,(t) being sampled at a 
frequency of M Hz, the normalized PSD Ey at a frequency f£ is defined as: 
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f 
J E(f )df 
En(f)=4 (2.3) 
[S(f) de 
0 
where E(f) is spectrum function for the signal I, (t), divided by the mean square deviation of I,. 
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Figure 7. PDFs of apparent source temperatures for smoldering fires. 
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Figure 8. The normalized PSDs of spectral radiation mtensities for the open fires. 


The spectrum function is the Fourier transform of the auto-correlation coefficient, R(T) 
and is given by: 


+co 
E(f)=2 [R(t)cos2xftdt (2.4) 


—0O 
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S(f) is the spectrum function for white noise, which when sampled at a frequency of M (and low 
pass filtered at M/2), has a uniform value 2/M between 0 and M/2, and is 0 at other frequencies. 


Physically, the normalized PSD at any frequency represents the ratio of the total energy of 
the signal below that frequency m comparison with a white noise source sampled with the same 
temporal resolution. The normalized PSDs shown im Fig. 8 imply that the total energy content of 
the spectral radiation itensity mcident on the detector below a frequency of 10 Hz for all the 
three open fires is approximately four time higher than that of white noise for direct radiation. 
For reflected radiation, smce the absolute values of the intensities are lower, the signal to noise 
ratio is lower. Therefore, the normalized PSD at 10 Hz is only about a factor of 3 higher than 
that of white noise. 


The normalized PSDs of spectral radiation intensities at 1000 nm for the two smoldermg 
fires are shown im Fig. 9. The mtensities obtamed from the direct view of the smoldermg wood 
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Figure 9. | The normalized PSDs of spectral radiation intensity for the smoldering fires. 


fire and from the reflected view of the cotton fire are barely above the noise level, and therefore 
could not be identified by the near infrared fire detector as valid fire signals. The major reason for 
the very poor performance with the smoldering wood fire is that the smolder surface is face down 
on the bumer and therefore very little of the smoldering surface is visible. The standard test fire 
(TE3) is specifically designed to test smoke detectors and not flame detectors. For flame 
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detection, an alternative standard test for smoldermg wood would be more appropriate. The 
smoldering cotton fire can be discrimmated utilizmg a A/D converter with a higher (16 bit) 


dynamic range. 


Based on the above results, the second criteria for the existence of the fire near the vicinity 
of the NIR fire detector was established. The second criteria is that the normalized PSD of the 
spectral radiation intensity at 10 Hz should be at least 1.5 time greater than that obtamed from 
white noise. This along with the first criteria (the 40% of the apparent source temperature should 
be between 600 and 2500 K) forms the basis of the fire detection algorithm of a NIR fire detector. 


The effect of different building materials on the bias m the temperature estimated by the 
NIR fire detector was checked by usmg a propane flame, and viewing the reflections using a 
window pane and a brick wall. The PDFs of temperature obtained from radiation reflected off the 
window pane and the brick wall is shown m Fig. 10. There is a potential for approximately a 200 
K upward bias m the temperatures estimated by a NIR fire detector. However, this still is well 
withm the temperature range specified by criteria 1, confirmimg the utility of the NIR fire detector 
for monitormg multiple room enclosures. 
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Figure 10. Effect of building materials on the apparent source temperatures. 
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The fmal part of the project converted the data acquisition and reduction program to 
enable the NIR fire detector to be used continuous monitoring. For continuously monitoring the 
environment for fire, data is collected for 25.6 seconds at 100 Hz providing 2560 intensity 
measurements at 900 and 1000 nm. From these tensity measurements a time series consisting of 
2560 apparent source temperatures were calculated and used to obtain the PDF and the average 
normalized PSD (from 20 digital transforms of 128 points). The data reduction program takes 
approximately 10 seconds on the laboratory computer. If the two criteria for the existence of a 
fire in the vicinity of the NIR fire detector was satisfied, an alarm was sounded. The three open 
test fires were lighted at different times during the data collection cycle, and the number of cycles 
required for the alarm to be activated was noted. This provides a rough estimate of the response 
time of the detector m units of cycle time and is shown in Table 1. 


Table 1: Response time for the NIR fire detector. 


The open heptane and foam fires were typically detected very quickly withm one or two 
data collection cycles (35 to 70 seconds) after ignition from both direct and reflected radiation. In 
addition, for almost all cases, the normalized PSD at 10 Hz was very distinct from white noise 
Discrimination by PDF is always successful as soon as the fire is started, smce the temperatures 
are very easily obtained. However, for the open wood fires, the fires had to reach a dimension of 
approximately 3 inches (for direct radiation) and 8 inches (for reflected radiation) before the 
normalized PSD was 1.5 times greater than that obtamed from white noise. In addition, no 
instances of false alarms were ever recorded during the entire operation over a two month period. 
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3. Numerical Evaluation 


3.1. DPF Formulation 


The primary objective of the numerical evaluation is to develop a methodology to 
calculate the imtensities incident on a fire detector from a source withm an enclosure. The 
anticipated benefit of the evaluation is to provide a design tool to help fire safety engineers 
determine the optimum location of a NIR fire detector for monitormg a multiple room enclosure. 
The spectral radiation mtensities incident on the fire detector withm an enclosure were obtamed 
using a ray tracing algorithm m conjunction with the DPF method. The basic numerical method 
was evaluated for two simple geometries (a cylindrical and a rectangular enclosure) durimg the 
present grant period. 


The first problem studied by the simulations was to determme the effect of the reflected 
radiation on the apparent source temperatures inferred by the NIR fire detector. For this purpose, 
a constant temperature source was placed at one end of a cylindrical enclosure as shown im Fig. 
11 (surface A,) and the detector was placed on the other end (surface Az). The cylindrical walls 


Source, A; A> (absorbing and reflecting surface) 


ate A3 


6 (Before 1“ reflection) 
6 (After 1* Reflection) 


Figure 11. Cylmdrical enclosure used to simulate the effect of reflections. 


of the enclosure was assumed to absorb and reflect (specularly and diffusely) the radiation 
incident on it. The spectral absorptivity and reflectivity of the cylindrical walls, in addition to the 
aspect (length to diameter) ratio of the enclosure, were the parameters of the simulation. 


The source emits photons im all angles and wavelengths. The number of photons emitted 
in a particular angle and wavelength depends on the directional and spectral emissivity of the 
source. The cylindrical surface of the enclosure, of length L and radius R, absorbs some of the 
photons incident on it and reflects the rest both specularly and diffusely. The trajectories of the 
photons emitted by the source and reaching the NIR fire detector were followed usmg the DPF 
method. The spectral radiation mtensity (i) incident on the detector at two wavelengths (A, and 
A2 ) were calculated and used to estimate the source temperature as: 
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where C3 is the second radiation constant and ¢, is the spectral emissivity. A temperature of 1500 
K and spectral emissivity values of 0.111 and 0.10 for the 900 and 1000 nm wavelengths were 
chosen for source. These values yield emission mtensities from the hot surface that are 
representative of the average intensities leaving luminous flames (Sivathanu et al, 1991). 


The surfaces and angles shown im Fig. 11 are discretized to obtam surface and angular 
ranges, similar to the method of Billings et al. (1991). The source surface A, is divided mto N 
surface ranges, which are of the shape of circular rmgs with equal width Ar each. The polar angle 
8, and the azimuthal angle o, are divided mto N equal angle ranges, AO and Ad respectively. 


The normalized probability density function (PDF) of a photon leavimg surface A; at a 
radius r, and polar and azimuthal angles 6 and 9 is given by (Sivathanu and Gore, 1994): 


f(r,n,0) =f(r)f(w)f@) (3.2) 
f(r) =2r/R? (3.3) 
f(u) = 2p (3.4) 
f(o)=1/2n (3.5) 


where is equal to cos 6 and R is the radius of the cylinder. Eqs. (3.2-3.5) represent Lambertian 
emission from the surface A;. Following Sivathanu and Gore (1993), the three PDFs associated 
with the emission of the photons from surface A, are discretized to obtain their respective DPFs 


P(r), P(), and P() as: 


P(r) = (33 Pi)31= LN (3.6) 

P(h) = (Hs Py), J=LN (3.7) 

P(o) = (O45 Pox )sk = 1,N (3.8) 
@; +A0/2 

where ieee heme, (3.9) 
@,-A®/2 


® represents one of the three variables r, 1, or @ and P,;, P,j and Py, are the probabilities that a 
photon is emitted from the surface range 1;, polar angular range 0; and azimuthal angular range Ox. 


The photon emitted from 1;, 0;, d, mtersects the cylmder at an x-location (x’): 
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where the superscript 1 represents the first mteraction of the photon with the cylindrical surface 
and J;, mj, and nj, are the direction cosmes given by: 


1; =Cos 0;;mj, =Sin 0;Cos 4505, =Sin 0;Sin $, (3.11) 


If x’. is greater than the length of the tube, the photon strikes the detector. Else, it intersects the 
wall where it undergoes a diffuse or specular reflection. The photons which leaves surface A; 
and intersects with the wall, will be characterized by a location given by the surface range x;' 
(which is on surface A), and new polar and azimuthal angles given by 0;', and ;', determined by 
whether it is specularly or diffusely reflected. This is identical to the two state parameters m the 
Markov chain method of Billings et al. (1991) and x;’, 0;', and @;,' can be considered to be the new 
three parameter state of the photon. The pomt of departure from the Markov cham method is 
that each of these states have a probability associated with them through the DPF calculations, 
rather than the probability bemg associated with the transition of a photon from the old state of r;, 
u;, and ¢, to the new state of x;’, Ly’, and o,.. 


The probabilities associated with the new location (x') of reflection treated as “Te- 
emission” of the photon are formed from all possible combmations of the old states and can be 
expressed as (Sivathanu and Gore, 1994): 


P(x!) = (xj';Py.)si= LN (3.12) 
N N N 1 

where = SFOS, Cbs Pri PujPorPa (3.13) 
i=1 j=l k=l 


where 8(x') = 1 if x; — Ax; /2<x! <x; + Ax; /2; 8x) = 0 else, and x! is obtained from 
Eq. (3.10). Physically, Eqs. (3.10-3.13) imply that the probability of a photon startmg from 
surface range x’; is the sum of all the products of the probabilities that the photon started from r;,, 
0;, and o, and mtersected the curved surface between x; - Ax;/2 and x; + Ax;/2. The probability of 
each photon is multiplied by the spectral reflectivity (p,) smce a certam percentage of these 
photons will be absorbed by the curved surface. Of the reflected a photons, a certain percentage 
is specularly reflected and the remainder is diffusely reflected. It is noted that for each value of i, j 


and k under the summation signs, the product 5(x!)P, PiyPok Px is identical to the transition 


probability of the photon goimg from states 1;, ;, and o, to the new state x’; after absorption by 
the surface. 
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It is slightly more difficult to calculate the probabilities for the angles of departure of the 
reflected photons, since these probabilities depend on the x-location. If the surface is totally 
diffuse (Stvathanu and Gore, 1994), then these probabilities are independent of the location x; and 
are identical to the probabilities of Lambertian emission. However, for the present problem, the 
surface reflects specularly and diffusely, and therefore, we need a jomt probability density function 
for the angles of emission and location of the photon after reflection. The method of calculating 
these jomt probability density functions are in principle identical to those given in Eqs. (3.10- 
3.13). For the polar angle 6’, and axial location x’, the joint discrete probability function of re- 
emission, is defined as (Sivathanu and Gore, 1993): 


eee tard 1 Jly.pl sips 
P(u »X )=((Hj>%j Ph ag )3j31=1LN (3.14) 


where Pi bor is the probability that a photon is emitted from the surface range x; in the polar angle 
Jol 


range 0; after 1 mteraction with the wall. For specular reflection the new polar angle of emission 
of the reflected photon which is incident at an angle of 0; is given by: 


(3.15) 


The probability associated with the photons re-emitted in the 6’; direction consists of both a 


specular component P Is Which is given by: 


1 jXj 
N N N 
PS => YY 8(x!)8@)PyPyPyP,-— = 3.18) 
MPM j=] j=l k=] (Pox + Par) 


where 8(0') = 1 if w; — Au; /2<Cos (6')<p; + Ap; /2, and 5(0') = 0 else. The 
superscript ‘s’ implies that only the spectral component has been used for calculatmg these 
probabilities. The last term in Eq. (3.16) represents the fraction of photons that undergoes 
specular reflection. The probability of the diffuse part is added to that of the specular component 
to obtain the joint discrete probability function of a photon emitted from the surface range x; m 
the polar angle range 0; after 1 interaction with the wall: 


pi, =pls = 2, ALS (Pap aa (3.17) 
b jXj HjXi a sk. 


The second term in Eq. (3.17) distributes the probabilities associated with the diffuse component 
to the different 0; direction based on Lambertian emission given by Eq. (3.4). 
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Finally, the joint DPF of axial location, polar and azimuthal angle P after 1 
kH jXj 


interaction with the wall is calculated in a similar manner. For specular reflection, the azimuthal 
angle of emission 6’ after one reflection, for a photon which was emitted from surface A; in the 
range Oy 1S: 


o' = by (3.18) 
and the associated probability is given by: 
N N N 
p 
pe =>) > Dy 8e')80 86 )PaP Pep, 
OEP Gy jal k=l (Psa + Pan) 


where 5(6') = 1 if 0, — Ady /2<o! <0, + Ad, /2, and 5(o') = 0 else. As before we 
add the probabilities associated with the diffuse component of reflection as: 


A N 
Pee = Pee A$ (pl yep ee (3.20) 
OE jx; OkHjXi 927 BX = pany KM FX 


At the end of the above procedure, Eq. (3.20) represents the probability associated the 
photon emission from surface A: (see Fig. 11) m the surface range x;, and angular range p;, and o, 
after one reflection. Summing up these probabilities for all the surface and angular ranges yields 
the total probability that a photon mitially emitted from surface A; will mteract with the surface 
A2. The remamimg probabilities represent the view factor (Sivathanu and Gore, 1994) from the 
hot end of the cylinder to the cold end. 


The next step im the process is to calculate the subsequent trajectories for the photons that 
have undergone one reflection with the wall. The radial location of the photons is no longer an 
mdependent variable, smce all the photons now start from the wall with rads R, however their 
axial location varies from 0 to L. The photons are started from all the N bins of x;’, p;', and ¢,’. 
The x coordinates of the second (and subsequent) mtersections are given by: 


2 
x? =x; +21; {Rmx / (m5. + Di )} (3.21) 


where the superscript ‘2’ is a counter mdicating the second interaction of the photon with the 
wall. The probabilities associated with each photon that undergoes this teraction are given by 


% . These probabilities are rebmned using the procedure explained in Egs. (3.13,3.17, and 

kB jX} 

3.20) to obtam p? . The solution procedure is carried forward to obtain p? : 
OkH jXi dich jXi 
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Bs , «ee. The process is termmated once more than 99.5% of the photons have been 
kH jXi 


absorbed by the cylmdrical surface or leave the enclosure through the right and left ends. 


The simulations to calculate the mtensity cident on the detector were carried out with 50 
ranges of polar and azimuthal angles and 50 surface ranges. For the Markov chain method, for 
even one interaction (e.g. the n-1 to the n™ interaction) the transition probability matrix would 
have 50x50x50 elements, with each element containing the probability of transition from x™’;, 6’, 
”', to x°;, 6°;, 6". These probabilities are difficult to compute when n is greater than 2 for such a 
large number of surface and angular ranges. In the DPF method, the probabilities associated with 
each interaction Eqs. (3.13, 3.17, and 3.20) were obtaimed by tracking a total of 125000 photon 
through just one interaction with the wall. It is noted that the recursive nature of the algorithm 
results m a computational advantage over explicit solutions, smce the probabilities are always 
rebinned after one mteraction with the wall, thereby restricting the total number of trajectory 
calculation to always 125000, despite the specular/diffuse nature of the mteraction and the three- 
dimensional trajectory. In addition, all possible photon trajectories (resolved to the surface and 
angular ranges) which have a probability of greater than 1 m 1000 are computed. This feature 
provides solutions with a high degree of accuracy. 


The disadvantage over the Markov chain method is that since the transition probabilities 
are not explicitly stored, any surface property changes will require the simulation to be repeated. 
However, if only the temperature of the source is changed, the probabilities can be explicitly used 
to find the intensity cident on the detector. 


The geometry used for the simulations m a rectangular enclosure is similar to that 
described for the axisymmetric enclosure. The rectangular enclosure of length L, width and 
height a, has an high temperature source (fire) at one end, and a detector at the other end as 
shown in Fig. 12. Part of the radiation emitted by the source is absorbed by the walls of the 
enclosure, and the remainder reflected both diffusely and specularly. 


Location of 
detector 


Location of 
fire 


Figure 12. Rectangular enclosure used to simulate the effect of reflections. 


To simulate a fire source at one end, the angular distribution of spectral radiation intensity 
leaving the fire is needed. This was obtained experimentally as described m the next section. 
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3.2. Description of Angular Measurements of Radiation Intensity 


The spectral radiation imtensities emanating from an heptane pool fire were measured at 
different polar angles using an optical fiber light guide as shown in Fig. 13. The objective of these 


Optical fiber 


Beamsplitter 


Filter (typ.) 
Heptane pool fire 


Ks 


Figure 13. Experimental arrangement used for the angular measurements. 


measurements was to obtain data to be used as mput for the numerical simulations mvolvimg a test 
fire within an arbitrarily shaped enclosure. The measurements were obtamed with a location on 
the axis at half the visible flame height as the reference pomt for the polar and azimuthal angles. 
The radiation emitted by the heptane pool fire was measured at 10 polar angles (8 from -60 to 
+60) using an optical fiber and two PMTs. The PMTs had optical filters centered at 900 and 
1000 nm (half bandwidth of 10 nm) m front of them The PMTs were calibrated usmg a 
blackbody and the jomt PDFs of spectral radiation mtensities at 900 and 1000 nm emanating from 
the fire for the different angles were obtamed. The heptane pool fire is axisymmetric. Therefore, 
the radiation is uniform im the azimuthal direction. 


A sample of the data obtamed from the heptane pool fire is shown m Fig. 14. The 
variation m the mean and RMS of spectral radiation mtensities with the polar angle is shown m 
the bottom panel of Fig. 14. As the polar angle either mcreases or decrease from horizontal, the 
mean spectral radiation mtensities increase. This is a direct effect of the mcrease m the path 
length with either an mcrease or decrease of the polar angle from zero. There is no variation with 
azimuthal angle for the axisymmetric (im the time averaged sense) pool fire. 


The variation of the mean apparent source temperatures with polar angle is shown im the 
top panel of Fig. 14. Despite the variation m mean spectral radiation intensities with the polar 
angle, the mean and RMS of apparent source temperatures are constant. This is to be expected 
since the temperatures estimated by the near mfrared fire detector are always biased towards the 
highest values along the path. 
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Figure 14. Variation of radiation properties with angle for an heptane pool fire. 


The DPF simulations described above were extended to accommodate the data obtaimed 
from these angular measurements of spectral radiation imtensities obtamed from the heptane pool 
fire. The heptane pool fire was assumed to be located at the center of one end of the cylindrical 
or rectangular enclosure. The fire detector was located at the other end of the enclosure. The 
aspect ratio of the enclosure was fixed at 3. The walls of the enclosure were assumed to have a 
coating with reflectivity of 0.8 at 1000 nm and 0.88 at 900 nm. The specularity of the reflectivity 
was set at 0.4. The spectral radiation mtensities obtamed from the angular measurements were 
input into the program as jomt probability density functions of spectral radiation mtensities which 
varied with the polar angle, 0. 
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3.3. Results of the Numerical Simulations 
The results of the simulations using a constant temperature source of 1500 K are 


presented first. The effect of the aspect ratio on the apparent source temperatures estimated by 
the NIR fire detector within a cylindrical enclosure is shown im Fig. 15. The walls of the 
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Figure 15. Estimated temperature variation with aspect ratio for a cylmdrical enclosure. 


enclosure were assumed to have an absorption coefficient which varied mversely with wavelength, 
typical of pigments used m commercial paints. When the aspect ratio is zero, the estimated 
temperature is 1500 K since only direct radiation is mcident on the detector. The absorptivity of 
the pamt was assumed to be 20% at 1000 nm.. As the aspect ratio of the cylmdrical enclosure 
mcreases, the radiation mcident on the detector consists of photons that have undergone multiple 
reflections with the walls of the enclosure. Therefore, the longer wavelength radiation is 
preferentially absorbed, leading to higher estimates for the apparent source temperatures. At an 
aspect ratio of 5, more that 70% of the photons have undergone at least 5 or more mteractions 
with the wall coating, and the measured source temperature is EE by approximately 500 K 
than the actual value. 


The effect of the specularity of reflectivity on the estimated temperatures is shown im Fig. 


16 for a cylindrical enclosure with an aspect ratio of 3. The effect of reflectivity depends on the 
number of imteractions the photons undergo with at the wall. An increase in the degree of 
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Figure 16. Estimated temperature variation with specularity for a cylindrical enclosure. 


specularity mcreases the number of photons that reach the detector after an interaction with the 
wall Therefore, the temperatures estimated by the detector are biased upwards since this larger 
fraction of the photons have undergone preferential longer wavelength absorption. When the 
specular reflection is unity, approximately 85% of the radiation incident on the detector is from 
photons that have undergone multiple reflections at the wall. 


The variation in the measured source temperatures with reflectivity of the cylmdrical 
surface is shown in Fig. 17. The ratio of specular to total reflectivity of the surface was set at 0.4. 
As the reflectivity of the cylindrical surface is decreased, a larger fraction of the photons is 
absorbed at the cylindrical surface and a higher percentage of the radiation mcident on the 
detector is from direct radiation, leading to a better estimate of the source temperature. 


The variation in estimated temperature with specularity of the reflectivity for a rectangular 
enclosure is shown in Fig. 18. The calculations were repeated usmg a Monte Carlo simulation 
with the same number of photons as in the DPF method . The accuracy of the DPF method over 
the Monte Carlo simulations for the same number photons is evident from Fig. 18. For highly 
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Figure 17. Estimated temperature variation with reflectivity for a cylindrical enclosure. 
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Figure 18. Estimated temperature variation with specularity for a rectangular enclosure. 
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specular surfaces, the effect of the reflected photons on the estimated source temperatures is not 
captured adequately by the Monte Carlo method. For highly specular surfaces, a larger fraction 
of the photons imcident on the detector has undergone multiple interactions with the wall, and 
adequate statistical sampling of these reflected photons (Sivathanu and Gore (1994))becomes an 
important issue. The Monte Carlo results would improve if a larger number of photons was 
initiated from the source. 


Unfortunately, there are no analytical solutions available that can be compared with the 
numerical results discussed above. However validation of the code was achieved by exammmg 
some specific conditions for which analytical solutions are available. When the reflectivity of the 
cylindrical walls was set to zero or 1 for both wavelengths, the temperature inferred by the 
simulation was 1500. In addition, when the reflectivity was set to zero, the normalized intensity 
incident on the detector was withm 0.01 % of the geometrical view factor. 


3.4 Results of Fire Detection Simulations 


Having validated the numerical code to a limited extent, the constant temperature source 
was replaced with the experimental data discussed in the previous section to simulate an heptane 
pool fire withm cylmdrical and rectangular enclosures. The PDFs of apparent source 
temperatures for a cylindrical enclosure with and without considering reflections from the walls of 
the enclosure are shown im Fig. 19. The walls of the enclosure were assumed to absorb 20% of 
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Figure 19. Effect of reflections on the apparent temperatures in a cylindrical enclosure. 
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the photons at 1000 nm, and 12% of the photons at 900 nm. The specularity of the reflectivity 
was also set at 0.4. A fraction of photons reaches the detector without undergomg any reflections 
from the wall. The apparent source temperatures estimated by the NIR fire detector from these 
photons are the same as those obtamed from direct viewing of the fire. A greater number of 
photons reaches the detector after interaction with the walls. This mteraction increases the 
apparent source temperatures estimated by the NIR fire detector smce longer wavelength photons 
are preferentially absorbed. Therefore, the PDF of apparent source temperatures obtained is a 
combination of direct and reflected photons, as indicated by the bimodal PDF shown mm Fig. 19. 


The effect of the wall reflections on the PDF of apparent source temperatures for the 
rectangular enclosure is shown im Fig. 20. Similar to the results obtamed within a cylindrical 
enclosure, the combmation of direct and reflected photons result m a bimodal PDF for the 
apparent source temperatures estimated by the NIR fire detector. 
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Figure 20. Effect of reflection on the apparent temperatures im a rectangular enclosure. 


Despite the higher values for the estimated source temperatures in both cases, the NIR fire 
detector could successfully discriminate the fires from background radiation., since most of the 
temperatures are still withm 800 to 2500 K. For a unit made of two rectangular shaped rooms, 
when the fire starts in one room, the above calculation will provide the jot PDF of intensities at 
the doorway to the other room, and also from the doorway to the any detector location. The 
present simulations helps to develop a design tool to optimally locate the fire detector in a 
multiple room enclosure. 
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4. Conclusions and Recommendations 


A new type of near infrared (NIR) fire detector was developed and evaluated during the 


two year period covered by this grant. The major conclusions of the present study are: 


(1) 


(2) 


(3) 
(4) 


(5) 


(1) 


(2) 
(3) 
(4) 


(5) 


(6) 


The normalized power spectral density of spectral radiation mtensity at a near infrared 
wavelength and the apparent source temperatures obtamed usmg two wavelength 
measurements are sufficient to determme the presence of open and smoldering fires. 

Open fires can be detected from direct as well as reflected radiation with very low false 
alarm rates using the NIR fire detector. Therefore, the NIR fire detector can be used for 
multiple room monitormg for open fires. 

Smoldermg fires can be detected only from direct radiation, smce the low temperatures 
associated with these fires resulted m very low radiation intensities at 900 and 1000 nm. 
The Discrete Probability Function (DPF) method m conjunction with a photon tracmg 
algorithm was extended to treat specular and diffuse reflections and absorption m a three 
dimensional enclosure. The DPF method is more accurate than the Monte Carlo method 
for determmimg the radiation mtensities mcident on the NIR fire detector after multiple 
reflections of the photons with the walls of the enclosure. 

The DPF method can be used to optimize the location of a NIR fire detector, used for 
monitormg a multiple room enclosure. 


The six specific recommendations for future work are: 


Characterize the spectral radiation mtensities at other wavelengths (particularly m the 
infrared) with a view of developing a low false alarm rate fire detector that can monitor 
both open and smoldering fires in multiple room enclosures. 

After the optimum wavelengths for fire detection have been determmed, continue the 
development of the algorithm to allow faster response times. 

Develop the DPF and other numerical technique for handlmg surface radiative heat 
transfer problems m other geometrical shapes. 

Study the various numerical techniques that can be used to Imk the results from simple 
geometrical shapes together so as to predict radiation behavior im multiple room 
enclosures. 

A Small Business Innovative Research (SBIR) program is necessary to develop a 
commercial prototype of this device. The requirements for commercialization include a 
fast programmable Digital Signal Processmg (DSP) Chip, CAD designed casimg for size 
optimization, and detector modification for current draw optimization. 

A second SBIR is needed for the imclusion of the DPF or Monte-Carlo based enclosure 
codes into architectural design software. 
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